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[57] ABSTRACT 

A method for determining the volume of compressible 
gas in a system including incompressible substances in a 
zero-gravity environment consisting of measuring the 
change in pressure (AP) for a known volume change 
rate (AV/At) in the polytropic region between isother- 
mal and adiabatic conditions. The measurements are 
utilized in an idealized formula for determining the 
change in isothermal pressure (AP^) for the gas. From 
the isothermal pressure change (AP^) the gas volume is 
obtained. The method is also applicable to determina- 
tion of gas volume by utilizing work (W) in the com- 
pression process. In a passive system, the relationship of 
specific densities can be obtained. 

7 Claims, 2 Drawing Sheets 
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VOLUMETRIC MEASUREMENT OF TANK 
VOLUME 

ORIGIN OF THE INVENTION 5 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-568 
(72 Stat. 435; 42 U.S.C. 2457). 

RELATED INVENTION 

This invention has subject matter related to the sub- 
ject matter disclosed in Ser. No. 217,725 filed on Jul. 11, 
1988 (MSC-21059-1). 

BACKGROUND OF THE INVENTION 

(1) Field of the Invention 

This invention relates to methods and apparatus for 2Q 
ascertaining the volume of liquid in a container or tank, 
and more particularly, to a system for gauging liquid 
volumetric quantities based upon the properties of a gas 
and is operable in a zero gravity environment. 

(2) Description of the Prior Art 25 

Prior Art Patents 

U.S. Pat Nos. 3,237,451, 3,413,847 and 3,769,834 all 
disclose methods to determine the volume of gases in 
tanks in a zero gravity situation. In U.S. Pat. No. 30 
3,237,451, an acoustic system is utilized for generating 
pressure changes which are related to volume measure- 
ments. In U.S. Pat. No. 3,413,847, density of a gas is 
measured and related to volume measurements. In U.S. 
Pat. No. 3,769,834, the volume of a human body is mea- 35 
sured by changes in pressure. 

U.S. Pat. No. 3,585,861 discloses a system for deter- 
mining volume of gases by using a reference gas pres- 
sure. U.S. Pat. No. 4,072,050 utilizes thermal energy 
changes for measuring volume. U.S. Pat. No. 4,384,925 40 
utilizes electrochemical sensing procedures for measur- 
ing volume. U.S. Pat. No. 4,726,216 relates to a calibra- 
tion circuit for detecting “HC” gases. U.S. Pat. No. 
4,781,061 relates temperature to volume measurements. 

BACKGROUND 45 

On the Earth, gauging systems for measuring liquid 
and gas quantities in a tank invariably rely on the 
Earth’s gravitational field to separate the gas and liquid 
phases in a consistent and predictable manner. This type 50 
of separation provides a flat, or predictable interface of 
the gas and liquid phases from which volumetric mea- 
surements can be made. Current state-of-the-art mea- 
surement instruments are designed to locate this separa- 
tion interface by some measured or inferred change in 55 
physical property. Knowledge of the storage tank ge- 
ometry is then used to calculate the respective volumes 
of gas and liquid. 

In a zero-gravity or weightless environment, how- 
ever, the gravitational or acceleration body forces are 60 
absent, which introduces a number of concerns due to 
the behavior of the individual phases, both dynamic and 
static, and the effect on contents gauging systems. 
These effects may include: 

1. no stable, known positioning or interface; configu- 65 
rations can take multiple geometries and this can 
result in multiple random distribution of interface 
locations; 
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2. for certain designs, wetting of sensor surfaces dis- 
ables measurement capabilities; 

3. loss of calibration. A calibration loss in a zero- 
gravity environment may make the system in-oper- 
ative since calibration is based upon a fixed, known 
geometry of the materials in the contained system; 

4. point sensor systems may not provide a true picture 
of the phase locations unless multiple sensors are 
used and statistical approaches are utilized. Contin- 
uous, total body systems require multiple detectors 
uniformly distributed over the container surface; 
and 

5. temperature and/or composition variations (strati- 
fication) can cause significant variations in mea- 
surement accuracy. 

Because of these effects, current methods for zero- 
gravity tank gauging are limited to single phase storage 
(either gas or supercritical) or rely on flow rate mea- 
surements into or out of the tank in the case of single 
phase storage, considerable weight is added to the 
structure since these systems must be designed for high 
pressure storage. 

Use of inlet/outlet flow gauging measurement may 
not provide a reliable (accurate) estimate of the tanks 
liquid contents due to the following restrictions: 

1. Two-phase flow. Most flow metering systems are 
designed to handle single phase systems. Two- 
phase metering designs are limited in application. 

2. Turn-down (range limitations). The turndown 
range of most flow gauging systems is at best 20 to 
1. Non-steady state applications of these instru- 
ments is significantly limited. 

3. System leaks. Storage tank leaks would go unde- 
tected since no direct means are available for con- 
tained liquid volume measurement. 

In the present invention, the contents in a two phase 
system including a compressible gas and an incompres- 
sible substance are measured independent of the shape 
of the container or tank and independent of tempera- 
ture. In another important aspect of the present inven- 
tion, the ratio of specific heat constants can be obtained 
in a relatively simple fashion. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a tank and measuring apparatus for 
practicing the invention; 

FIG. 2 is a plot of change in pressure versus rate of 
volume change; 

FIG. 3 is a schematic representation of an analog 
computer system for equation solving; and 

FIG. 4 is a schematic representation of a digital CPU 
unit for equation solving. 

DESCRIPTION OF THE PRESENT INVENTION 

The present invention may best be understood by 
reference to underlying principles utilized in the inven- 
tion. 

In a storage tank having a predetermined volume and 
containing a gas and liquid, the gas is compressed and 
the compressibility properties of the gas are used to 
provide a measurement of the liquid volume of the tank 
by virtue of the following relationship. 

V Uq= V tank ~ V gas (0 

where V// ? is the volume of the liquid, V tank is the 
known volume of the tank, when empty and V gas is the 
volume of gas. 
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The gas volume V^is determined by measuring the 
total pressure in the tank and the pressure required to 
change the volume of the gas in the tank by a small 
known or predetermined amount in a defined period of 
time. The measurements are utilized in a gas compress- 5 
ibility relationship to determine the gas volume and 
hence the liquid volume relative to the known tank 
volume. 

In the present invention, a formulation was devel- 
oped to utilize the measurements of a change in pressure 10 
(AP) where the volume of the gas is changed by a fixed 
value (AV). From the pressure (AP) measurements, the 
gas volume (V gaj ) can be calculated by use of the fol- 
lowing relationship. 




—NP&V 

A/» 


Where “N” is a compressibility constant which is 
dependent upon the amount of heat generated during 20 
compression; P is static pressure in a tank; AP is the 
measured change in pressure; and AV is the predeter- 
mined change in volume. The value of “N” varies from 
“1” for an isothermal condition and a maximum of 
about “1.6” for an adiabatic condition for some gases. 25 
Most gases of interest to the present invention fall 
within the range of 1 to 1.6 for “N” 


where P is pressure, V is volume and N is a value depen- 
dent upon the amount of heat transferred during the 
compression process. Equation 3 is the equation of 
change for a gas compression process between the limits 
of N equal to a value of one (1) for isothermal condi- 
tions to a value about 1.67 for adiabatic conditions. The 
upper limit of N at adiabatic conditions is equal to ratio 
of specific heats (K) which is the ratio of Cp/Cy where 
Cp is the gas heat capacity for constant pressure and Cy 
is the gas heat capacity for constant volume. The ratio 
of Cp/Cyis sometimes referred to as the ratio of specific 
heats. This ratio ranges between approximately +1.1 
and +1.67. It has been determined by analytical and 
experimental means that the value of “N” between the 
isothermal region and the adiabatic region is a function 
of the rate of heat transferred from the compressed/ex- 
panded gas to its surroundings, such as the tank wall. 
The rate of heat transferred is dependent upon the time 
rate o volume change (or volume displacement rate). 
Thus, the value for "N” is bounded between two ther- 
modynamic regions where the lower limit of “N” is an 
isothermal compression region (constant temperature) 
where the value of “N” is about 1. The upper limit of 
“N” is the adiabatic compression region (no heat trans- 
fer) where the value of “N” is equal to the ratio of 


Referring first to FIG. 1 for orientation relative to 
physical apparatus, a tank 10 of a known volume V 
contains a compressible fluid and incompressible sub- 30 
stances. A motor 11 drives a volume displacement de- 
vice 12 which may be a bellows, diaphragm, piston etc. 
that is driven in a known linear or oscillating manner to 
produce a known volume change (AV) in the tank. The 
resulting pressure change (AP) is a function of the com- 35 
pressibility factor “N” of the gas. The change in volume 
can be measured by a AV volume measuring device 13. 

A total pressure gauge device 14 measures the mean 
pressure in the tank and is coupled to the tank 10 but 
isolated from instantaneous pressure changes by a snub- 40 
ber 15, if required. A AP pressure measuring device 16 
measures the change in pressure. The motor displace- 
ment 11,12 produces a constant value of change in vol- 
ume AV in a time period At where the time At can be 
varied for the known change in volume AV. 45 

With a given gas, the gas can be compressed in an 
isothermal compression region or condition (the com- 
pression is accomplished by a change in volume over a 
very long time period) where the temperature remains 
constant. The given gas also can be compressed in an 50 
adiabatic compression region or condition (the com- 
pression is accomplished by a change in volume over 
very short time period) where there is no measureable 
heat transfer during the change in volume. 

In the present invention a gas is compressed in a tank 55 
with a known change of volume (AV) for a defined 
change of time (At) between the critical points for an 
isothermal condition and an adiabatic condition so that 
the value of the isothermal pressure change AP^can be 
determined. The gas volume in the tank can be deter- 60 
mined once the value for AP;*, is determined. As set 
forth in Equation 2, with a known change in volume 
AV, for a given change in time At and a given pressure 
(P) and a measured change in differential pressure (AP) 
at isothermal conditions where the value of “N” is one, 65 
the volume of the gas (Vgas) can be determined. Equa- 
tion 2 was derived by taking the derivative of 


specific heats K. The ratio of specific heats K is the ratio 
of constant pressure specific heat to constant volume 
specific heat. The ratio K is a function of gas compres- 
sion and operating conditions, but as noted above, is 
usually in the range of 1.1 to 1.6. 

The foregoing described relationship is illustrated in 
FIG. 2. In FIG. 2, the isothermal region 20 extends to a 
polytropic region 21 to the adiabatic region 22. In the 
present invention, the characteristics of the polytropic 
region 21 are utilized to determine the volume of gas in 
a tank dependent upon the type of gas, the volume 
displacement ratio, and the heat transfer characteristics. 

While precise characteristics of a given gas between 
the isothermal region 20 and the adiabatic region 22 can 
vary, the precise shape of the curve in the polytropic 
region 21 generally will have the form illustrated and 
will vary in scalings with repect to the AV/At scale. 
Data was obtained and correlated by statistical methods 
to determine the shape or for the curve 25 for a given 
gas and a given tank. As shown in FIG. 2, at low vol- 
ume displacement rates (region 28), environmental ef- 
fects dominate the shape of the curve and the value of 
AP will tend to converge to a limiting value in the 
isothermal region 20 as the volume displacement rate is 
increased. Volume displace rate is synonymous with a 
frequency of volume change. At the limiting value of 
the isothermal region 20, the value of “N” is equal to 1 . 
In the polytrophic region 21, the response of AP, as a 
function of an increase in the volume displacement rate, 
is a monotonic functional response. In the adiabatic 
region 22, “N” is equal to "K”. The rate of volume 
displacement increases to a point 29 where error in 
measurements occur due to slow response time of sen- 
sors, vapor/liquid dynamic interference, accoustical 
noise generated by the volume displacement device and 
so forth. Measurements above the point 29 may be unre- 
liable because of equipment noise. 

The characteristic of the curve 25 in the polytropic 
region 21 is generally defined by the following equation 


AP/=AP (M (Ae -«//+!), 


(4) 
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Where AP / is the measured differential pressure for 
the i th measurement, AP^) is the differential pressure 
for isothermal compression (a constant), fi as the fre- 
quency (AV/At), and A and a are constants. From a 5 
numerical solution of this equation for A and a, the 
value for AP^) can be predicted. Solution for the vol- 
ume is then given by the following relation: 




-PAK 


(5) 


where N=l. 

To obtain the value of AP^ the preferred method is 
to operate the motor 11 at each of three discrete fre- 
quencies (AV) and for each frequency to determine the 
change in pressure AP. The change in pressure AP^ can 
be found by use of Equation 4. With the values of AV, 

At and AP, these measurements can be used to deter- 20 
mine the gas volume in the following described manner. 

Equation 4 above, can be restated as step function as 
follows: 

AP f *APuo(Ae-"+l) (6) 25 

where 

AP / is the initial change in pressure, 

APzjo is the stabilized differential pressure under iso- 
thermal conditions, 30 

A is a constant which is equal to (Cp/Cv— 1) where 
Cp is the specific heat of a gas at a constant pres- 
sure and Cv is the specific heat of a gas at a con- 
stant volume, 

e~ a/ is a time dependent variable which goes from 1 35 
to 0 as time (t) increases. 

A rapid volume displacement causes the immediate 
pressure pulse which decays according to Equation 6. A 
Computer or other suitable instrumentation performs a 
numerical analysis for solution of the constants P/^ A 40 
and a in the equations. From AP/s 0 the volume of gas 
can be determined from Equation 5. 

In further detail of the above, the time dependent 
variable, when redefined as a time period T of sinusoi- 
dal change of pressure, is 45 

(7) 
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The term P ^ 0 is a computed variable which is deter- 
mined as follows: Using the notation that 

P itt/Vf-X ( 11 ) 

three equations can be written representing each of the 
three frequency measurement periods, as follows: 


= X^e-^ + 1) 
= X(Ae-W + 1 ) 
= XlAe-W + 1 ) 


( 12 ) 


With the time period T^less than time period T 2 and 
the time period T 2 is less than T 3 , three equations and 
three unknowns A, fi and P/^can be solved as follows: 


— - x^XAe-en 
V\ 


— - X = XAe~P n 

n 


— - X = XAe-W 
v 3 


03) 


The terms A and are eliminated by mathematical 
simplification and by substitution, the following rela- 
tionship is obtained. 


0)-* [(f)- 
(t)-* [($•)- 


( 73 — T2)/(TZ — Tl) 


(14) 


The solution of Equation 14 for the value of X can be 
made by analog or digital computer techniques during 
each measurement cycle to give a continuously derived 
value of X that satisfies Equation 14. 

Since 


When the value of AP/ is changed to a cyclic value, 
then the value of AP/ becomes a fluctuating pressure 50 
around the pressure value P which is the average pres- 
sure. The pressure change AP/ has time period T/. To 
use this value in computation, it is divided by the total 
volume change and is a non-cyclic variable. Using the 
following notations for ease of explanation, where the 55 
notations are 


— A V /ae V,AP, = p/ 


( 8 ) 


the non-cyclic variable is 


60 


P//v/ 


(9) 


When corrected for phase shift, this gives 


EL 

V/ 


Piso 


{Ac-W 


1) 


65 

00 ) 


P(iso) A Piso —p (15) 

= vi * ~"£7“ = iwr 

and with a known total system pressure, the volume of 
gas in the tank can be calculated by 

V(gas)=P/X (16) 

Based upon actual measurement values, the differ- 
ence in the calculated gas volume determined by the 
present invention to the actual gas volume in a tank was 
3.71 percent. 

In utilizing this method it should be noted that the gas 
volume determination is independent of the properties 
of the gas in the tank. 

The ratio value of specific heats C p/C yean be easily 
obtained with the present invention. The value of 
Qp/Cy is equal to tank volume times the change in 
pressure times the tank pressure divided by the change 
in volume. 
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Referring now to FIG. 3, an analog computer system 
is schematically illustrated. Pressure readings of pres- 
sure are supplied by a pressure sensor 14 in electrical 
signal form to a computation network 30. The differen- 
tial pressure sensor 16 supplied differential pressure 5 
measurements in electrical signal form to a correlation 
amplifier 31. The differential volume device 13 is oper- 
ated by a linear actuator 32 which also drives a linear 
encoder 33 to provide an electrical signal which is pro- 
cessed by a counter 34 and a digital to analog converter 10 
35 to a correlation amplifier 36. The correlation ampli- 
fier 36 supplies the volume measurements for three 
different times. Separate oscillators 40-42 drive the 
linear actuator 32 at selected constant volume changes 
and rates. The output of the oscillators 40-42 also is I 5 
supplied to the correlation amplifier 36 and to a period 
to voltage conversion circuit. The output of the period 
to conversion network is supplied to the computation 
network 30. 

The volume displacement device for these systems is 20 
on the order of 1/1000 or less of the total tank volume. 
Experimental studies have shown that this volume ratio 
can be 1/100,000 without significant loss in accuracy. 
The minimum volume change is determined by the 
sensitivity of the differential pressure and volume dis- 25 
placement sensors. The rate of volume change is typi- 
cally in the cycles per second range. But this frequency 
may vary over several orders of magnitude as depen- 
dent on the sophistication of the measurement system. 
The volume displacement device will nominally be of 30 
an expanding metal bellows design, since the reliability 
(life) of these systems are for 1,000,000 cycles or more. 

Typical differential gauges used in experimental stud- 
ies have measurement capabilities that range from 0 to 
0.1 in H 2 O (total span), response times are listed as 0.5 35 
milliseconds. The differential pressure gauge, both in 
minimum range and in response time, sets the major 
design parameters for this device. 

As shown in FIG. 4, a CPU 40 can be programmed to 
receive the AP, AV/AT and pressure measurements for 40 
three frequencies of operation and digitally process the 
data for determining the volume of gas. 

Another approach to determining the volume of gas 
in a tank is to obtain values for pressure (P) and volume 
(V) at two different frequencies. If y is a known value of 45 
the ratio of specific heats C/>/CKand the term 


(7-0* A 


( 17 ) 


is substituted for A in the first two equation sets of 50 
Equation 13 then the following equation sets result. 


£L 

n 


( 18 ) 


X(y - 1) 


— — X 

n 

- i) 


: e -/3n 




55 


The term e~0 is eliminated by combining the two 
equation sets of Equation 18 which give 


60 


/ \ ✓ x 71/77 

Or)-* ft-*) 

{X(y - 1)> “ K*(y ~ 1 ) ) 


( 19 ) 
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The solution of Equation 19 for the valued X is simi- 
lar to the solution of Equation 14 above. 

In another aspect of the invention, a work meter 45 as 
shown in FIG. 1 derives an electrical signal representa- 
tive of work during the compression process and is used 
to sense the work input “W” during the compression 
process so that “N” can be determined from the equa- 
tion, 


^ML, where ^, (20) 

After determining “N’\ the value of “N” is substi- 
tuted into Equation 2 for determining the gas volume. 
As before, the formulation solution can be a computer 
using programming in a conventional manner. 

The system of the present invention further monitors 
the density change (or temperature or pressure) for a 
given volume change so that the gas volume can be 
determined. This is basically a passive tank gauging 
measurement since no complex moving parts are re- 
quired. The driving force that causes the volume 
change is due to variations in the tank temperature due 
to external environmental temperature changes. In 
other words, a temperature change external to the tank 
produces a pressure rise or fall in the tank. 

The relationship of the volume of a simple compress- 
ible substance for most gases is 


V(gas) 


— AV 

1 - (V2/V1) 


(21) 


For most gases, PVT equations of state can be used 
for solution of specific volume which gives 


^ {gas ) 


— AV 


1 - 


1 ± A T/T \ 
1 4- A P/P ) 


(22) 


Thus with the determination of temperature changes 
and pressure changes over a long period of time, the gas 
volume can be determined from the relation set forth in 
Equation 22. 

It will be apparent to those skilled in the art that 
various changes may be made in the invention without 
departing from the spirit and scope thereof and there- 
fore the invention is not limited by that which is en- 
closed in the drawings and specifications but only as 
indicated in the appended claims. 

We claim: 

1, A method for determining the volume of a con- 
tainer containing at least one compressible gas and at 
least one incompressible substance independent of the 
shape of the container and independent of temperature 
by the steps of: 

producing a change in volume (AV) and pressure 
change (AP) in said container over a known period 
of time (At) and determining the change pressure 
(AP) in said container while said gas is in a poly- 
tropic region between a compressibility factor (N) 
equal to one (1) for an isothermal region and a 
compressibility factor (N) equal to K, the adiabatic 
coefficient for an adiabatic region, 

establishing a predetermined relationship of the 
change in pressure AP as a function of the rate of 
change of volume (AV/At) for said polytropic 
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region between said isothermal region and said AP/=AP&XAe a '+i) 


adiabatic region, 

determining from said predetermined relationship, 
the value of the change in pressure (AP^) for iso- 
thermal conditions, and 5 

determining the volume of the compressible gas by 
measuring the median gas pressure (P) and utilizing 
the measured gas pressure P, the change in volume 
AV and the determined AP /*> in the formula 

10 

v -PAY 

V{gas) “ A/W) * 

2. The method as set forth in claim 1 wherein said 
predetermined relationship is defined by the formula 15 

AP/=AP (iw) (Ae-“^+l) t 

where P/is the measured differential pressure for the i th 
measurement, P^is the differential pressure for isother- 20 
mal conditions (a constant), fi is the frequency (Av/At), 
and A and a are constants. 

3. The method as set forth in claim 2 wherein the AP/ 

determination is made for at least three distinct values of 
AV/At in the polytropic region. 25 

4 . The method as set forth in claim 3 wherein the 
change in volume(AV) is held constant and the value of 
At is varied to obtain the three different values of 
AV/At. 

5. The method as set forth in claim 4 wherein the 20 
determination of P, AP/ and At for said three different 
values are solved by simultaneous equations to obtain 
the volume of gas. 

6. A method for determining the volume of a con- 35 
tainer containing at least one compressible gas and at ** 
least one incompressible substance independent of the 
shape of the container and independent of temperature 
by the steps of: 

producing a change in pressure (AP) value and a ^ 
change in volume (AV) in said container over each 
of three different time periods in a range of values 
for AV/At between an isothermal condition and an 
adiabatic condition, 

determining the change of pressure AP^ for isother- 45 
mal conditions from the relationship of 


where 

At is the time period 
AP/ is the initial change in pressure, 

AP/so is the stabilized differential pressure under 
isothermal conditions, 

A is a constant which is equal to (Cp/Cv— 1) 
where Cp is the specific heat of a gas at a con- 
stant pressure and Cv is the specific heat of a gas 
at a constant volume, 

e~ a 'is a time dependent variable which goes from 
1 to 0 as time (t) increases, 
for each of the three time periods, 
determining the volume of gas by the use of the deter- 
mined value of AP/s 0 in the relationship of 

r/ _ —PAY 

V(gas) ~ a/wT • 

7. A method for determining the volume of a con- 
tainer containing at least one compressible gas and at 
least one incompressible substance independent of the 
shape of the container and independent of temperature 
by the steps of: 

dynamically producing a known change in volume 
(AV) and pressure change (AP) in said container 
and determining the adiabatic pressure, determin- 
ing the subsequent isothermal pressure in said con- 
tainer after said gas reaches an isothermal pressure, 
and determining the time period to change from 
adiabatic pressure to isothermal pressure, 
then determining the volume of gas from the relation- 
ship 

AP/=AP/jo(Ae~ af + 1) 
where 

AP/ is the initial change in pressure, 

A Pjso is the stabilized differential pressure under 
isothermal conditions, 

A is a constant which is equal to (Cp/Cv — 1) where Cp 
is the specific heat of a gas at a constant pressure and Cv 
is the specific heat of a gas at a constant volume, 
e~ a/ is a time dependent variable which goes from 1 
to 0 as time (t) increases. 

* * * * * 
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